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' Abstract 

o . 

■ The phenomenon of macroscopic homogenization is illustrated with a simple example of 
! diffusion. We examine the conditions under which a d-dimensional simple random walk in a 
' symmetric random media converges to a Brownian motion. For d = 1, both the macroscopic 

■ homogeneity condition and the diffusion coefficient can be read from an explicit expression 
for the Green's function. Except for this case, the two available formulas for the effective 



diffusion matrix k do not explicit show how macroscopic homogenization takes place. Using 
an electrostatic analogy due to Anshelevich, Khanin and Sinai we discuss upper and 

lower bounds on the diffusion coefficient k for d > 1. 



1 Introduction 

X 

d [ The long time behavior of random walks on a random environment is reviewed. We focus mainly 
on the following question: 

What are the conditions under which a properly scaled random walk on a non-homogeneous 
medium converges to a Brownian motion. 

This and related phenomena are usually named macroscopic homogenization (and the envi- 
ronmental conditions are called macroscopic homogeneity conditions) because such system looks 
homogeneous at macroscopic scales. The discussion will be restricted to simple random walks with 
non-vanishing transition probabilities (or rates) {w^xy)} such that 

W{xy) = W(^y^) (1.1) 

holds for all nearest neighbor sites (xy) of a d -dimensional lattice Z'^. The so-called symmetric 
medium has been considered by several authors (see e.g. jABSUt lAKSt lAVt IMFGW| IKu| IPVj and 
references therein). 

*e-niail: marchett@ime.usp.br 
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Except for a basic lemma, the general scheme of our presentation will be dimensional indepen- 
dent. However, the one-dimensional problem plays a central role in this work since, in this case, 
the macroscopic homogeneity condition can be read from an explicit formula. 

The d = 1 case has been mostly investigated. The first mathematical results |KKS[ ISo[ ISi] were 
concerned with asymmetric random walks with transition probabilities Wx,x+i = 1 —Wx,x-i, a; G Z, 
being independent and identically distributed (i.i.d.) random variables (note that Wx,x+i 7^ uix+i,x)- 
The trajectories {X{t),t > 0|X(0) = Xq} of asymmetric random walks were shown to behave 
very anomalously. Symmetric random walks began to be discussed in a series of papers (see e.g. 
|ABU1 iBSWl lABSUn in connection with the problem of disordered chains of harmonic oscillators 
(see |LMj for an introduction and a selection of reprints) and other problems in physics. Using 
Dyson's integral equation 0, these authors derived the following asymptotic behavior for the 
trajectories: if 

E^;,i.+i<oo, (1.2) 
where E denotes the expectation value with respect to the i.i.d. random variables {wx^x+i}, then 

as t — > 00, and the convergence of X(t)/\/t to the Gaussian random variable with zero mean and 
variance 2 {E,w~^^j^^ ^ is implied. In addition, if condition ()1.2p is violated, then EX^()f:) grows as 
with an exponent 5 <1 depending on the divergence of the distribution at Wx^x+i = 0. 

A mathematical proof of convergence to Brownian motion for c? = 1 symmetric random walks 
was given by Anshelevich and Vologodski |AVj . For d >2 there are at least two different proofs and 
both require macroscopic homogeneity conditions more stringent than ()1.2|1 . Anshelevich, Khanin 
and Sinai [AKS^ proved the result by developing an expansion for the expected value of the inverse 
of a non-homogeneous discrete Laplacian. Kiinnemann |Kuj has proven this result by extending 
Papanicolaou- Varadhan's approach |PVj . Whether p.2|) is a necessary (and sufficient) macroscopic 
homogeneity condition for d >2, remains, to our knowledge, an open problem. 

The present paper is inspired by the work of Anshelevich, Khanin and Sinai. We use the logic 
of this proof in order to simplify the Anshelevich and Vologodski's proofs for d = 1. Our proof, in 
particular, eliminates the technical hypothesis of Wx,x+i to be strictly positive and illustrates with 
textbook's mathematical methods the macroscopic homogeneity condition ()1.2|) . 

A Brownian motion is described by the diffusion equation. The random environment induces 
an effective diffusive matrix k, whose elements are given by limj__>oo E {Xi{t)Xj{t)) /t, i,j = 1, . . . ,d. 
For the one-dimensional problem, the reciprocal of this constant is given by the macroscopic homo- 
geneity condition: k,~^ = Kw~l._^^i. For d >2, the two available formulas for the effective diffusion 
matrix (see |AK^ IKuj ) do not explicitly show how macroscopic homogenization takes place and 
this makes it difficult to obtain estimates. In this review (see also (AKSij ) the upper bound 

K < EW{xy) 

will be shown to hold if |l — W(^xy)/'^uj(xy)\ < 5 < 1/2. We also discuss how a lower bound can be 
obtained using the electrostatic equivalence of the diffusion problem as formulated in |AKSj . 

The outline of the present work is as follows. In Section |21 we formulate the problem and state 
our results. The proofs will be given in Section El under the assumption that the eigenvalues and 
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eigenvectors of the semi-group generator of the process converge to the eigenvalues and eigenfunc- 
tions of a Laplacian. The eigenvalue problem is a consequence of our basic lemma (Lemma I3.6p 
which will be proven in Section |3]for d = Ihj Green's function method. The spectral perturbation 
theory will be presented in Section El Finally, the diffusion coefficient will be examined in Section 

El 



2 Statement of Results 

Let B denote the set of bonds of the regular lattice Z"' and let w = {wh}beM be an assignment of 
positive numbers. Each component Wb represents the transition rate of a random walk to go from 
the site x to y along the bond b = (xy). The assignment w defines an symmetric environment on 
Z*^ if the transition rates satisfy W{xy} = wi^yx) ■ 

Given an environment w and a finite set A C Z"^, a continuous time random walk on A, with 
absorbing boundary condition, is a Markov process {XA,io(t), t > 0} with differential transition 
matrix PVa = Wp^iw) defined by 



W{xy) Uy 



\x-y\^l 



{Aa,wu)^ + u^ ^ W(^y), (2.3) 



\x-y\^l 



where 



XJ 7 



\x-y\^l 



for all a; G A and u such that = if y G Z'^ \ A. 
Note that — Aa,^ is a positive matrix, 

{u, (-Aa,^)u) = ^ ^ W(^^y) {uy - Un,f > , (2.4) 
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x,yezd. 

\x-y\ = l 

and, if Wb = w for all 6 G B, Aa,«j = w Aa where Aa is the finite difference Laplacian with 0- 
Dirichlet boundary condition on A. From here on, A is taken to be the hypercube centered at origin 
with size |A| = (2L — l)'^, L G N: Al := {x = {xi, . . . , Xd) G Z"' : supj \xi\ < L}, and all quantities 
depending on A will be indexed by L instead. 

The probability distribution of {XL^uj{t), t > 0} is governed by the semi-group generated 
by Al^w If XL^w{t) denotes the position of a random walk at time t, then 

P{XL,Ut) = x\Xl,M = 0) = (e*^^-'Mo)^ , (2.5) 

where (uq)^ = 5o,x- 

The semi-group is the solution of the initial value problem in M^, 



■t 

du 



Ai,^u = 0, (2.6) 
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with initial condition u{0) = Uq. 

The solution of ()2.6p exists for all times t > and all sizes L < oo but may depend on the 
realization of w and on the initial value. We present the sufficient conditions on the environment w 
by which the solution of ()2.6j) . under suitable scaling of time and space, converges almost everywhere 
in w to the fundamental solution of the heat equation, 

|-#„,^0. (2.7) 

with u{t, &D) = 0. Here, ()2.7j) is defined in the domain t > and ^ G P := (— 1, 1)'^ with boundary 
dV = : supj l^il = 1}, and = d'^{n) is given by 

d' = d-Kd=J2'^^j^- (2-8) 

i,j=l 

The heat kernel Tt{ri, Q = e^'^^^^r], Q when defined in x M'^ gives rise to a Wiener process (or 
Brownian motion) {B{t), t > 0} with covariance E,B{t)B{s) = 4fi;min(s,t) (see e.g. Simon jS]). In 
view of the boundary condition, Tt yields a Brownian motion {Bolt), t > 0} on the domain T> with 
absorbing frontiers. The probability density of Bq{s + t) — Bq{s) to be equal to C, can be obtained 
by solving equations (j2.7|) by Fourier method 

d 

m,0)= e-(-'/^)"'*n^^(™^^'/2) ' (2-9) 

where = n - nn = Ylt j=i^i'^ij^i sc(njx) stands for either cos (n^x) or sin (riix) depending on 
whether rii is an odd or even number. 

This discussion suggests the following definition. 

Definition 2.1 The random walk in a random environment w is said to converge to a Brownian 
motion if there exist a constant k = k{w), the diffusion coefficient, such that 



hm (e^^*^--«o) = (e^^'^'^Vo) (0 



m 

uniformly int > 0, ^ &V and /io on the space of finite measures M. with support in T> . For r G M'^, 
[r] G if" and has as components the integer part of the components of r; and Uq G M"^^ is a vector 
given by 

Axi+l)/L i-{xa+l)/L 

(wo). = L'' ■■■ d/xo(C) • 

Jxi/L Jx^/L 

It is important to note that, by the dominant convergence theorem, this definition implies the 
convergence in distribution of the random walk process {{1 / L)X L^w{L'^t) , t > 0} to the Brownian 
motion {i?o(t), t > 0} as it is known in Probability Theory (weak convergence of their distributions): 

X^it) = — „,(L^t) — > Boit) in distribution if, for any collection < Si < ■ ■ • < s„ of positive 
Ij 

numbers and any collection /i, ...,/„ of bounded and continuous functions in n G N, we have 

E^o/i(^l(Si)) ■ • • fniXLiSn)) ^ E^Ji(5o(Sl)) " " " /n(i?o(Sn)) , (2.10) 
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as N ^ oo, where E^^ means the expectation of the process starting with the measure /xq- Note also 
that Xilt) has been scaled as in the central limit theorem: ^^(t) is a sum of about independent 
increments^ divided up by L. 

Theorem 2.2 (Anshelevich and Vologodski )AVj ) If d = 1 and the environment w is a sta- 
tionary process such that the partial sums 



converge as x ^ oo to k ^ , < k < oo, almost everywhere in w, then a random walk in a random 
environment w converges, for almost every w, to a Brownian motion with diffusion coefficient n. 

Theorem 2.3 (Anshelevich, Khanin and Sinai jAKS| ) For any d > 1 and 5 < 1/2, let w = 
{wb}beM be independent and identically distributed random variables such that 



with w =Ewb. Then a random walk in a random environment w converges, for almost every w, to 
a Brownian motion with diffusion coefficient matrix k, = K,{d,6,w). 

Remark 2.4 1. Theorem \2/J was proven in lAVf under the condition \2.11\) with Wx,x+i > 0. 
The positive condition has been eliminated in our proof. Condition \2.11]) is met if w = {wb\ 
are i.i d. random variables with < 'Kw^'^ < oo. Finiteness of first negative moment seems to 
be, according to arguments presented in lABSOf (see also IFINf ). a sufficient and necessary 
homogeneity condition since, otherwise, the walk would spend a extremely large time between 
jumps leading the process to be sub-diffusive. 

2. Whether the homogeneity condition < Eiw^^ < oo is also sufficient for d > 1 is, to our 
knowledge, an open problem. It would already be an important achievement to prove Theorem 
\2.!^ for any distribution whose support is M. It is unfortunate that both proofs (see lAKSl U^jlf ) 
require as a technical step Wb to be bounded away from and oo . 

3. Theorem \2.'J\ holds also if the transition probabilities {wxy} do no vanish for y — x belonging 
to a finite set iV~ that is able to generate Z*^ by translations. Under the same sort of condition 
\2.12^) . \ AK!:^ have shown convergence to Brownian motion satisfying the diffusion equation 



4. Theorem \2.'^ can also be extended to Brownian motion on if one combines the result with 
both absorbing and periodic boundary conditions on &D (see details in lAKSf ). 

simple random walk with continuum time jumps according to a Poisson process on time with rate 2c? and 
there will be 2dLF't jumps in average after a time V^t. With the random environment, the Poisson process has a site 
dependent rate given by Y.y:\x~y\=i^{xv)- 






(2.12) 
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Theorem 2.5 Under the conditions of Theorem \2. and a conjecture formulated in \6.21^ , there 
exist a finite constant C = C{d), such that the bounds 

1-Q<-^<1 (2.13) 

hold with 1 being the d x d identity matrix and g a positive matrix such that 

450(1 + C) 
^ - 1 - 4(5 (1 + C) ' 

in the sense of quadratic forms. 



3 Basic Lemma 

The proof of Theorems 12 . 21 and 12 . 31 presented in this section are based on the uniform convergence 
of the eigenvalues and eigenvectors of A^^^ to the eigenvalues and eigenfunctions of the Laplacian 
operator on V. The uniform convergence follows from a classical result in perturbation theory 
which says the following. 

If y4i, ^2, . . . , v4„, . . . is a sequence of bounded operators in a Hilbert space TC which converge to 
A in the operator norm, then all isolated pieces of their spectrum and their respective projections 
converge uniformly, as n ^ oo, to those of A. 

Because and d"^ are unbounded operators we consider their inverse instead. To formulate 
the results of this section, we need some definitions. 

Let ii be an isometry of the vector space into the piecewise constant functions in the vector 
space Lq(T'), of square-integrable functions / on "D = (—1, 1)'^ with f{dV) = 0: 

given by 

'^11^:'^ (3.14, 

with [x] as in Definition 12.11 

The adjoint operator : Ll(V) — ^ R^, 

dxdfix), (3.15) 

/L Jxa/L 



Jxi/L J X, 



is defined by the equation 

{f ,iLu) = {i[f ,u), 
with the inner product in and Ll{V) given, respectively, by 
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and 



{f,g)= [ f{x)g{x)dx. (3.17) 
We shall denote by 

E2'■.= L-'{i,All,^i) (3.18) 

the scaled inverse of A/,^^. The operator kernel of is a step function with step-width 1/L 
which, as we shall see in the next lemma, approximates the kernel (9^)~^(r, s) in the operator norm 
induced by Lg-norm: 

\\A\\:= sup P/II2, (3.19) 

/:||/Il2 = l 

where := (/,/). 

Lemma 3.6 (Basic Lemma) Under the conditions of Theorems \2. 'A and \2. there exist a finite 
number Lq = Lq{w) such that, {S^"^;L > Lq] is a sequence of hounded self-adjoint operator in 
Ll{V) which converges 



-1 (92)-i||^0, (3.20) 



II^L 

as L —>■ 00, in the operator norm topology. 

It thus follows from perturbation theory (see e.g. [F]): 

Corollary 3.7 // A is an isolated eigenvalue of (i9^)~^ and E the orthogonal projection in its 
eigenspace £, one can find a subspace £^ C Lq{T)), invariant under S^^, and a corresponding 
orthogonal projection such that 

\\E^-E\\ — ^0 and \\E^ {El^ - Xl) E^\\ — >0, (3.21) 

as L ^ 00 

Lemma f3. 61 will be proven for d = 1 random walks in Section 0] This lemma reduces the Brow- 
nian motion limit problem to the convergence of the inverse matrix A^^^ to the inverse Laplacian 
((9^)~^ in the Lg-operator norm topology. Corollarv 13. 71 will be proven in Sectional 

Proof of Theorems \2.S\ and assuming Corollary \3. 1\ (As in Appendix 3 of |AKSj ) Let 
: ]R_ — > M be given by 

^{X) = e*/^ (3.22) 
and note that ip is uniformly continuous at A = with ip{0) = 0. We have 

Tt = ^iid')-') and T,^ = e*^- = ^ (S^^) . 

In view of Definition 12. II and the isometry ii, Theorems 12.21 and 12.31 can be restated as^ 

sup sup I [{Tt - r,^)/io] (0 1 ^ , (3.23) 



The isometry has been introduced to bring aU operators to the ScLiiie Hilbert sp&-CG Z/q(2?). Not6, however, 
I' and T,^ 



that ^ and remain finite rank operators 
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as L — s> oo. We shall prove an equivalent statement: 



lim \\Tt 

L— >oo 



lim M{d'r') ~ ^ {^t) 11=0. 



L-^oo 



The inverse Laplacian (9^) ^ on P with -Dirichlet boundary condition is a compact operator 
with spectral decomposition given by (recall equation ()2.9p ) 



(3.24) 



-1 



where 



TT/ 



r/2)2 ^ 

eigenvalues and associate eigenf unctions of (9^)~^ and 



4/(7r2n2) and e„(0 = Y[sc[{7i/2)niii\, n e N^, are 



i=l 



Because is an accumulation point, we introduce an integer cut-off N < oo and let 

(^n)^^ = ^ An En . 



(3.25) 



(3.26) 



n|<JV 



We have 



|n|>Ar 



(3.27) 



\n\>N 

which can be made as small as we wish by letting — > oo. More generally, the uniform continuity 
of V9 at means the following: given e > and a non-positive bounded operators A, we can find 
6 > such that if ||A|| < 5 we have ||y9(y4)|| < e/3. We shall use this fact often in the sequel. 

From Corollarv 13. 7[ there exist a projector 



E 



L,N 



|n|<JV 



where E!^ is, analogously to En, the projector on the invariant subspace: '^i^Sn = £^ 



^L.N 



L 



N) 



Writing 
(3.28) 

(3.29) 



and using the fact that E ' is an orthogonal projector, we have 

(here ip{A~^) is defined by its power series I + tA + t^A^/2 + ■■■). 

We now show that — can be made smaller then 5 = 5{e) if L and are chosen 

sufficient large. By Lemma EIEl there exist Li > Lq, Li = Li{6), such that 



< 



6 



(3.30) 



■^)-^-S-V||<^ (3.32) 



for all L> Li. From (IT^ and (IT^ . there exist A^i = Ni{5) such that 

\\idV-{dl)-'\\<^-, (3.31) 
if > A?"!. By Lemma ESI there exist L2 > Lq, L2 = L2{S), such that 

3 

holds for all L > L2 with fixed. 

If L > max(Li, L2) and > A^i, equations ^^-^^ yield 

W^l' - ^l!n\\ < W^l' - (dY'W + Wid'r' - id%)-'\\ + 11(4)-' - ^l!n\\ < 5, (3.33) 
which implies, due the continuity of if and the orthogonal relation ()3.29p . 

M^L)-vi^L,N)\\<'-. (3.34) 

By uniform continuity of ip and ()3.32j) . we also have 

M{dV)-^{{dl)-')\\<'-. (3.35) 

In addition, using the spectral decomposition of (<9^)~^ and S^^jy, and taking into account 

||A„,-E„ — Xn^nW ^ l^n — ^nlW^nW + ||-E„ — (3.36) 

and Lemma (3.61 we can find L3 > Lq, L3 = L^{e, N), such that 

Midi)-') - ^{^-^],)\\ < I . (3.37) 

Now, let L > max(Li,L2,L3). It then follows from dSISl, (HH^ and (pn?7|l 

+ ymr') - ^(s^vijii + m^,') - ^(s^y ii < e , 

which implies strong convergence of the semi-group and completes the proof of Theorem 12.31 □ 

Remark 3.8 The introduction of the cut-off' N in \3.2(^) is necessary even for homogeneous envi- 
ronment. In this case, the eigenvalues and eigenvectors of L^^ Aj^^ can be computed explicitly: 

(d d 
AL^ sin^ ^n, j and e^(x) = JJ (^^^^ ' (3-38) 

i=l / i=l 

With n e A* := {1, 2, . . . (2L - 1)}'^ (recall sc (uix) stands for cos (uix) or sin (n^x), depending on 
whether rii is an odd or even number). Note that |A„ — A^|, with A„ given by \3.24\ ), may not be 
small if \n\ = 0{L). We always pick N large but fixed and let L 00 in order ^3.37\ ) to be true. 
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4 Proof of Lemma 13.61 for d = 1 



In this section we prove Lemma 13.61 for d = 1. We consider the second-order Sturm-Liouville 
difference equation 

^L,wU = f 

with u{dV) = 0, and use the method of Green to calculate the matrix elements of A^^^. This gives, 
in view of equation ()3.18|) . an explicit formula for the operator kernel (r, s). 

The procedure starts by looking for two linear independent solutions of the homogeneous equa- 
tion 

(Al,^ u)^ = ?i'(x-l,x>(«x'-l - U^) - W{x,x+l){Ux - Ux+l) = , (4.1) 

with X G {— L + 1, . . . , L — 1} and = ul = 0. Without loss of generality, we set tf(L-i,L) = 

W{-L-L+l) = /«• 

Proposition 4.1 Let C,l ^ M^^^"*^ be a vector valued function of the environment w given by 

X 

y=-L+l 

for all X G {—L + 1, . . . , L — 1}, where 

L 

y=-L+l 

Then Ui = C,l o-nd U2 = 1 — C,l ore two linear independent solutions of ^-l]) . 

Proof. Ul = is a solution of ()4.H1 by simple verification and the same can be said of U2 = 1 — Cl- 
For this, note that 

W{x-i,x)i'^^L)x = Vl (4.3) 

holds uniformly in x, where {'Vu)^ = — u^-i- It thus remains to verify that they are linear 
independent. 

Let W = W{ui, U2] x) be the "Wronskian" of the two solutions Ui and U2 given by the following 
determinant 

iUl)^ {U2)x 
W(x-1,x>(VMi)x W(a'-l,a;>(VM2)x 

It follows that two solutions Ui and U2 are linear independent if W{ui,U2;x) ^ for all x G 
{— L, . . . , L}. Plugging Ul and U2 into ()4.4|) we have, in view of ()4.3|) . 



W 



(4.4) 



W = -r]L [iU)x + 1 - (a).] = -VL , (4.5) 
which concludes the proof of the proposition. □ 
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The inverse matrix can be calculated by the so called Green's function method (see e.g. 



x,y 



-Vl 

[1 - (q).] iU)y 

-Vl 



if X < y 



if X > y . 



(4.6) 



To see this is true, we note (^Aj^^^)^^ is the 2;-component of a vector for each y fixed. So, by 
definition 



holds for all X ^ y. For x = y we have 

/T-7,- \ {^LjX fTlt- \ ^ \sl)x 

= -W(x,x+1>(V4l)x+1 W(x-1,x)(V4l)x'- 



(a)x + (1 - (a).) = 1 



by (j4.3j) . verifying the assertion. 



We are now ready to write the operator kernel of . In view of 

{z^g,Ai^f) = J 9)xA,^y{z^ f)y 



ix+l)/L 
x,yeA -^x/L 
1 /.I 



1 r(^+'-)/-^ riy+'^)/L 

- V L / drg{r)A^^yL / ds^( 



(4.7) 



1 ^-1 



dr I dsg{r) [L AyLr\\Ls\) f{s) = {g,iAi^ f), 



valid for any (2L - 1) x (2L - 1) matrix A and / G Lg(P), and definitions (j3T8|l and (jTUjl . we 
have 



[LrULs] ' 

for any — 1 < r, s < 1 . 
If new variables 

Cl.= 2^l-1 

are introduced into equation ()4.(i|l . the operator kernel ()4.8|1 can be written as 

-1 



(4.^ 



ir,s) = — — {1 - |(CL)[Lr] - {Cl)[Ls] \ - iCL)[Lr] (Cl)[Ls]} 



4Lr] 



(4.9) 



in view of the fact that (Cl)x is a monotone increasing function of x. 
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By Schwarz inequality the operator norm ()3.19|) is bounded by the L^-norm of the operator 
kernel, the Hilbert-Schmidt norm IHi^'Hp := [ dr [ ds\K{r,s)\'^. Since the functions in the 
Hilbert space has compact support, we have 



1 j-i 



ll^ll < llli^lll < 4 sup |i^(r,s)| 

-l<r,s<l 

and — — > is implied by the pointwise convergence S^^(r, s) — > {d^)~^{r,s) of 

the operator kernel. We shall see that the latter convergence sense is consequence of the following 
result. 

Proposition 4.2 Given e > and w satisfying the hypothesis of Theorem \2//A Then, there exist 
an integer number Lq = Lo{e,w) such that 

\{CL)[Lr]-r\<e (4.10) 

holds for all L > Lq and — 1 < r < 1. 

Proof. Under the hypothesis of Theorem 12.21 the strong law of large numbers holds and 

2Lr]L-^n, (4.11) 
as L — * oo, for almost every w (see eq. ()4.2j) ). Analogously, since [Lr]/L — > r as L oo, we have 

/ > N \Lr] + L 1 i-4 -, 

(C^)[-^']=2^^^-^^-TZhTT S ^' (4-12) 

for each r G (—1, 1) and this gives ()4.1()j] . □ 

The Green's function method can also be used to compute the integral kernel of (9^)~^ as an 
operator in the Hilbert space Lq(D). The two linear independent solutions of the homogeneous 
equation 

/€--4 = 0, -l<r<l, (4.13) 

with boundary condition /(—I) = /(I) = are /i = 1 + r and /2 = 1 — r. Replacing Ui(2) and 
W(x-i,x)C^Ui(^2)) by /i(2) and i^{dfi(2)/dr) in ()4.4|) . gives W = —2k. Substituting these into ()4.6|) 
following the simplification of (I4.9|l . yields 

{d^)-^{r,s) = —{l-\r - s\-rs} (4.14) 

Note that \{d'^)~^{r, s)\ < 1/{2k) and, in view of (03) and Proposition |Ol 

\^l\r,s)\<C (4.15) 
holds uniformly in r, s G (— 1, 1) for all L > Lq. 
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Now, let pL(r) := {i(L){r)-r and := 2LriL^l^/K. Then, if L > Lq, in view of (gH), 
and Proposition 14. 2| we have 



^(r,s)-{dY\r,s)\ < 



2Lr]i 



- 1 



S^^(r,s)-(a^)-^(r,.) , (4.16) 



where 



S^(r, s)-(92) ^(r, s) = — |pL(r)pL(s) + rpL(s) + spL(r) + Ipl(s) - pL(r) 



2k 
2e 

< \ 

K 2k 



(4.17) 



uniformly in r, s G (—1, 1). When combined with 1)4.111) and ()4.15)1 this proves Lemma ITHl for d = 1. 

□ 



5 Perturbation of Spectra 

Proof of Corollary \3. 7\ This proof can be found in Appendix B of |AKSj and is essentially based 
on the perturbation theory of Hermitian bounded operators developed by Friedrichs in iFj. Since 
it can be described shortly, we repeat the proof's derivation for completeness. Our derivation, 
however, is more close to in the sense that we perturb an interval of the spectrum. When the 
interval contains one single eigenvalue this reduces to the derivation of [AKS|. The generalization 
is however essential in dealing with intervals containing accumulation point of the spectrum. This 
situation has to be considered in order to show that the spectrum projection in such intervals 
remains orthogonal when the perturbation is turned on. 

We now introduce some notation. Let Jq G M be an isolated closed interval of the spectrum 
(t((9~^) of defined with Dirichlet boundary condition on T> = (—1, 1). There exist < S < oo 
and an interval / C /q such that / fl <j{d~'^) = (j{d~^) fl Jq and 

dist (Jo, R\I)>S. 

Let Sq denote the eigenspace spanje^ : A„ G Jq} G LI(V) associated with Jq and Eq the spectral 
projection onto Sq. Let £i/l denote the subspace of the Hilbert space TC := Ll(V) invariant under 
the action of S/, and Ei/l the projection (not necessarily orthogonal) onto Si/l. 

The projection Ei/^ is defined by the following set of equations: 

{I-Ei/L)El'EyL = (5.18) 

(i.e. Si/L is an invariant subspace) and the two conditions 

EyLEo = EyL and EoEyL = Eo. (5.19) 

Note that, under ()5.19j) -Ei/l is a projector 

^l/L = i^l/L Eo)Ei/L = Eiil{Eq Ei/l) = Ei/L Eq = Ei/L , 



13 



which is consistent with Eq in the sense that hmi_+oo -E-i/l = -£"0- 

We shall prove that, provided Vl := — d^"^ is bounded, i?i/L depends analytically on 1/L 
and £i/L tends to £^0 as 1/L ^ 0. The proof of this statement uses equation (jS.lSj) to write an 
integral equation. For simplicity, we shall drop the index L of the quantities H^^, Vl, £i/l and 



Ei/L- 



Our stating point begin with equation 

{1-E)d-^E = d-^E-Ed-^E 



E - Ed' 



2 



(5.20) 



which comes from the following facts. 

The operator 9"^ commutes with the spectral projector Eq. Using this and equations ()5.19|) . 
we have 

Ed-^ E = E Eod'^ E = Ed'^ EqE = E O'^Eq = E E^d'^ = E d'^ , 

and this implies the second line of ()5.20|) . The commutation relation = allows us to 

replace E in the equation ()5.2U|) hy Q := E — Eq 

{1 - E)d-^ E = d-^ Q - Q d-\ (5.21) 

Combining ()5.18|1 with ()5.2H1 and using = 9"^ + V, gives 

d-'Q = Qd-'-iI~E)VE 

= d-^Q-{I-Eo-Q)V{Eo + Q). ^^-^^^ 

Since the interval /q is isolate from the rest of the spectrum, d~'^ is an invertible bounded 
operator in the subspace (J — Eq) H. We can solve the left hand side of ()5.22|1 for Q by defining 



Xf:- 



d'f iifeil-Eo)n 

otherwise . 



Note that X d'^ = d'^ X = I - Eq and ||X|| < 6-\ 
Equation ()5.22|) can thus be written as 

Q = g{Q), (5.23) 

where 

giQ)=Xd~''Q-X{I + Eo-Q)ViQ- Eo) . 



Proposition 5.3 The sequence Qn, n = 0, 1, . . of projectors defined by 

Qn = 9iQn-l) , (5.24) 

with initial condition Qo = satisfies the conditions i5.iy\} and converges, Q = lim^^ocQn; to the 
unique solution of equation 45. 2!^} . 
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Proof. We have ||Qn|| < <? < 1 for all n G N provided q is chosen small enough and L is taken so 
large that if < q then 

||^(Q)1| < ||X|| ll^o^^'ll q + (1 + qf\\VL\\ < q. (5.25) 
Note that the smallness of g depends on the smallness of V. Since 

\\X\\\\Eod-^\\<-^^:=a<l, 

equation ()5.25|) holds provided 

11-11 (-) 

Now, for fixed value of q, it can be shown (see ref. [F| for details) 

\9{Q)-9{Q')\<0\Q-Q'\ 

also holds with 6 < 1 and this implies Proposition 15.31 by the Banach fixed point theorem. □ 

We have proven the existence of a unique projector Ei/^ such that ||Qi/l|| = II-E'i/l — Eq\\ < q. 
Since q can be made arbitrarily small by taking L sufficiently large so that ()5.26|) holds, we have 
limj;^_,oo \\Ei/L — EqW = and Si/^ — > Sq as L ^ oo. 

To complete the proof of Theorem 13.61 we need to find an orthogonal projector E^ onto = 
Si/L in order to get ()3.21|) . This is achieved by setting 

= E,/LiEl^E,/L)-'Elr. (5-27) 
and noting that the inverse operator (E'^E)'^ exist because \\Q\\ < q implies 

\\E,f\\<\\Ef\\+q\\f\\ = \\Ef\\+q\\Eofh 
for any f E H such that E^f = f. As a consequence \\E f\\ > (1 — q)\\f\\ and 

{f,E^Ef)>{l-qnf\\. 

One can show, in addition, that [S^^, E^] = for all L. Therefore, for any A G /q H o"((9~^), 

\\E^ {E'^ - XI) E^W = \\ {E-^ - XI) E^W 

< \\ (d-^ - XI) {Eo + {E^^ - Eo)) \\ + \\V E'^W (5.28) 

< II (9-2 - XI) II + \\d-^\\ \\E^ - EoW + \\Vl\\ . 

Since the right hand side goes to zero as L — > oo this concludes the proof of CoroUarv 13.71 □ 
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6 Diffusion Coefficient 



This section is devoted to the proof of Theorem 12.51 The diffusion coefficient will be estimated 
throughout an expansion for the expectation of the inverse matrix, E(Al^^)~^, with w satisfying 
the macroscopic homogeneity condition ()2.12j) . This is justified in ref. [AKS^ in view of the fact 
that {Al^w)~^, when properly scaled, converge to its expectation for almost all environment w. 
Thus, the formula (E(Al^^)^^) ^ ~ Al,k, is expect to hold in the limit as L ^ oo. We will see that 
very important cancellations take place by inverting the series expansion of K{Al,w)^^ ■ 

A simple algebraic manipulation shows 

{-^L,^f" {-^L,^f" = T^f^ , (6.1) 

where 

Dl,^ := (-A^,^)-^/' {Al,^ - Al,w) i-AL,wy'^' (6.2) 

is a well defined matrix since, in view of (j2.4j] and (j2.12|l . — Al,«j is positive and the square root of 
—Al;w can be taken. 

Choosing W = Kwh and use ()2.3p to write Aj^_^ = WAi where A^ is the finite difference 
Laplacian with 0-Dirichlet boundary condition on A, equation (|6.2j) can be written as 



Dl,^ = (-Ai)-^/2 Ai,, (-Ai)-^/2 , (6.3) 
where a = {ab} given by at = Wb/w — 1, are i.i.d. random variables with mean Ea^ = 0, such that 

\ab\ <5<^ (6.4) 

holds in view of (jZT^ . 

Equation ()6.1|) suggests us the use of Neumann series to develop a formal expansion of (A/, „,)~^ 
in power of -Dl.ui due to the small parameter S. The remaining of this section is devoted to the 

pointwise convergence of the matrix element of [E (/ — Dl,w) ^] 
Using ()6.3|) . we have 

= / + E (-^^)"'' [(-^-)" ^^-] (-^-)"'' • (6-5) 

' n>l 



To write (j6.5|) in a more convenient form, let Vl : — > M'^^ be the finite difference operator: 

(VlM)(xJ/> = -(VLM)(y^) = ai^^y) [Uy " M^) . 

where the sign ci^xy) = YliiVi ~ ^i) — =^1, according to whether (xy) is positively (= 1) or negatively 
(= — 1) oriented. Vl maps a 0-form u into a 1-forms VlU- Let V2 : ^^'^ — ^ be its adjoint 
{uj, Vl u) = (V2 1^, u), i.e. the finite divergent operator which maps a 1-form uj into a 0-form V^ciJ 
given by 

(V2 uj)x = E ^{y^) ' 

y:\x-y\=l 
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and let : M®-^ — > M^'- be the multiplication operator by a: {MaUj)h := abUJ^. 
With these notations, we have 

Al,„ = V2M,Vl, (6.6) 

and its bilinear form reads 



recovering expression ()2.4|) for the quadratic form. 
Define 

$:= Vl(-Al)-WI, (6.7) 

and note that, since (— Ai)~|^ is the Coulomb potential between two unit charges located at x and 
y, ^h,b' is the dipole interaction potential between two unit dipoles located at b and b'. Note that $ 
maps 1-form into 1-form. 

In view of fj6.6p and ()6.7|) . equation ()6.5p can be rewritten as 



^ = ^ + EE"rW^r, (6. 



where, for T = (6i, 62, • • • , &n), 

n 

ar = ]^a6fe, (6.9) 

fc=i 

and ^ 

(i;, VTrw) = ^^61 $62,63 ■ ■ ■ '^b„.ubn ^br. , (6.10) 

with and being 1-forms given by V (— A^)""^^^ v and V (— A^)""'^''^ u, respectively. 

Concern the convergence, as A Z"^, of a generic term of the expansion ()6.8p . the following 
remark is now in order. 



Remark 6.1 The asymptotic behavior of the dipole potential ^b,b' for L » dist (&,&') >> 1 can 
be estimated by its spectrum decomposition ^ , 



where 




(6.11) 



k=l 



^Since we have not rescaled the space A C Z'', it is convenient to introduce a base {e^}„gA» , A* = {1, . . . , 2L— 1}'*, 
normaUzed with respect to the scalar product ((m,i')) := X^^eA "a: ^'a: = L'^{u^v): = e^/v'T^. The spectrum 
resolution of the identity is written in terms of this base. 
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and as in 113. 3^) . If we take b = {xx^'^^) and b' = {yy^-'^) , where z^^"^ is a nearest site of z 
whose components are given by z^^^ = Zi + 6k/, and make a change of variables, = {tx /2L)ni, 
i = 1, . . . ,d, we have 

lim = / rfV f Vsin2(<^,/2) I V*Vmcos(<^fc(a;fc-yfc)) , (6.12) 



where f{z) = f{z^'^^) — f(z) is the difference operator in the k-th direction. The \x — y\ >> 
1 behavior of ^h,b' is given by restricting the integral id.l!^} around a e -neighborhood of with 
e\x — y\ = 0(1); 

^b,b' ~ / d'^ifi fj-f± tan {^i[xi - yi)) tan (v5j(xj - yj)) TT cos (<^fc(xfc - yk)) 

V-^r J\^\<e ^J^ 

dh^-^ (6.13) 



{2'KY\x-y\'' J\t\<o(i) 



[dist (6, b')Y' 

As a consequence, ^b,b' is not summable in absolute value, 

\^b,b'\ ~ logL, 

and the uniform convergence with respect to A of the r-summation in ()6.8j) requires cancellations 
due to the dipole orientations (see ref. jPPNM]). 

We shall exhibit in the following another kind of cancellation due to the inversion of the expected 
value of 

Inverting the expectation of ()6.8|) gives 

[E{i~DL,u,yy = i-eL, (6.14) 

where 

0^ = E(-1)'^' E E Ear, •••Ear.W^r, (6.15) 

k>l ni,...,nfcr=(ri,...,rfe) 

and r G X ■ ■ ■ X B^N with rii > 1. Note that Wr = Wr,--- Wr,. 

To see how the log-divergent terms in ()6.14|) cancel out, it is convenient to use graph-theoretical 
language. A graph G consists of two sets {V,E): V = {vi, . . . ,Vs} is the vertex set and E = 
{ei, . . . ,es'} the connecting set of edges. To each edge e its assigned an ordered pair of vertices 
{vv') (its extremities) which are called adjacent if f 7^ v'; otherwise e is said to be a "loop". To 
the problem at our hand, we shall identify the bonds {bi, . . . , 6„} as a the vertex set of a graph G 
whose connectivity is determined by the presence of interactions ^b,b'- 

Two graphs G and G' are isomorphic (denoted G ~ G') if there is a one-to-one correspondence 
between their elements which preserves the incidence relation. A path F on G is an ordered sequence 
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{f jg, Cj^, . . . , ei„, f j,^} of alternately vertices and edges of G such that Cj^. = {yi^_^Vi^^ holds for each 
k] the edges {cj^, . . . , Cj^} are the steps of the path and the vertices {w^q, . . . , t>j^} are the points 
visited by the path. V may be identified with one of these ordered sets since it can be uniquely 
determined by each of them. Two vertices v,v' &V may be connected by more than one path. A 
graph G is said to be connected if any two vertices f , v' can be joined by at least one path F on 
G. The components of a non-connected graph are its maximum connected subgraphs. Given two 
vertices the disconnecting set of edges is a set whose removal from the graph G destroys all 
paths between v,v'. A cut-set is a minimal set of edges the removal of which from a connected 
graph G causes it to fall into two components Gi,G2- 

Turning back to equation (j6.14p . one may interpret V = {bi, . . . , 6„} as a set of vertices visited 
by a path. In view of the fact that ai, has zero mean, we have 

Ear = ar = (6.16) 

if there exist at least one bond bi which are not repeated in the list F = {61, . . . , 

The condition ()6.16|) says that the path F must visit each vertex at least twice otherwise 
its contribution to (j6.14p vanishes. The set of distinct bonds V = and edges E = 

{(6162) , • • • , (bn-ibn)} form a connected graph G with even valency V(6) > 4 for each vertex b E G. 
Graphs with this property will be called admissible graphs. Note that each path F yields only one 
graph G but there are possibly many n-step paths covering each edge (bi-ibi) of G exactly once 
which starts at bi and ends at 6„. If we denote by [F]g the set of all paths F satisfying these 
conditions for a given admissible graph G, we have 



Proposition 6.2 Equation AS. 15]) can be written as 

Ql = Y1 Y1 "^gWg, (6.17) 



n>l G:\E\^n-l 
admissible 



with 

Ag:= E (-ir'«r. ■■■«r., (6.18) 

re[r]Gr=(ri,...,r,) 

where we sum over all sizes n G N, n > 1, all admissible graphs G of size \E\ = n, over all paths 
F in [F]g and over all decompositions of V into s,s > 1, successive paths (Fi,...,Fs), each of 
which capable of generating admissible graphs Gi. Here, with the notation of \6.10\) and footnote 
in Remark \6.1l 

{{v,Wgu))=A, I n "^^A '"^r., (6.19) 
\(bb')eEiG) ) 

for n> 1 with ^biM = 1 for n = 1 (the case that G is the trivial graph ({61}, 0) ). 



We shall in the sequel state two lemmas and prove Theorem 12.51 under an extra assumption. 



Lemma 6.3 If G is an admissible graph with at least one cut-set contained one edge (i.e. G falls 
into two components by cutting a single edge), then Aq = 0. 
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Lemma 6.4 There exist a constant Cig] < oo, depending on the equivalence class [G] of isomorphic 
graphs G, such that 

E 

bi ,bn fixed 

holds uniformly in L for all admissible graph G with \E\ = n and cut-sets with no less than two 
elements. 



n 

{bb')eE(G') 



< 



C[G] 



dist(6i, bn 



i2d 



(6.20) 



Remark 6.5 The proof of Lemmas \6.'^ and "KJ are essentially given in lAKS^ (see Assertions I and 
II of Section 4)- Note that our estimate i6.2(J\) have not included the logarithmic corrections which 
appears in that reference. To get rid of these one has to control the loop subgraphs of G carefully 
as it is done in the ref. VFPNM l. The uniform upper bound \6.2U\} results from the hypothesis that 
G remains connected by cutting one single edge. 



Graphs with single edge cut-sets do not contribute to ()6.17j) due to the following cancellation 
in Lemma (6. H[ 

Proof of Lemma 16.31 Let {pi &j+i) be the only edge of a cut-set and let F = (Fi, . . . , F^) be a 
decomposition of a path in G. Either both fej and belongs to some Fj or they belong to two 
successive ones. We call the latter decomposition type A and the former type B. It turns out 
that there is an one-to-one correspondence between type A and type B decompositions differing 
only by the splitting of F^ into two elements F^^"* and F^-^'*. Lemma (6.31 follows from the fact that 
the contribution to ()6.18|) of a pair of decompositions established by this correspondence have the 
same absolute value with opposite signals. Note ar = a {i)a {2) if the edge (6j&i+i) bridges the 

two paths F^^'' and 



□ 



In view of Proposition 16.21 and Lemmas 16.31 and 16. 4t equation ()6.15|) can be estimated as 



b,b' 

where 



Kl^Y. E -^IGiqol, (6.22) 



n>l [G]:\E\^u-l 
admissible 



with the sum running over the equivalence classes [G] of isomorphic admissible graphs G of size 
\E\= n and G\g] as in Lemma f6.4[ Note that Ag = ^[g]- 

Now we show that, if one uses, as in refs. |AKSj and |PPNMj . the upper bound 

C[G] < (6.23) 

for some geometric constant C < oo where r = \V\ is the number of vertices in G, the equation 
(I6.22|) cannot be bounded uniformly in L. Taking into account property (|6.4I) . 

|ari ■ ■ ■ arj < 5" (6.24) 
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holds uniformly in T and equation ()6.22p can be bounded by 



/r, _ 1 \ min{n,|A|) 

^^^E \ E ^(n,r)C^ (6.25) 

n>l ^ ^ r=l 

Here, we have identified each path F = {61, . . . , bn} in a given graph G = {V, E) of size \E\ = n — 1 
with a partition P = (Pi,...,Pr) of {1? 2, . . . , ra} into r = \V\ pairwise disjoint subsets. This 
association is one-to-one since F is an ordered set of elements. Note that each component Pj 
corresponds to one vertex bi^ = = ■ ■ ■ = bi^ of G visited p = \Pj\ times by the path F. One 
can thus replace the sum over all equivalent classes of graphs [G] and over all paths F in [G] by 
the sum over all partitions P. The factor n(n, r) counts the number of partitions of {1, 2, ... , n} 
into r subsets. The bound (j6.25j) disregards the fact that G is an admissible graph. Also, the 
consistency of each decomposition F = (Fi, . . . , Fg) into paths Fj's which gives rise to admissible 
graphs has been not considered. The binomial factor (^"^^) < 4" counts the decomposition of F 
with n steps into any number of paths with the number of steps < n (the cardinality of the set 
{1 < "^1 < "^2 < ■ ■ ■ in-i < n — 1}). In addition, for the upper limit in the second sum we note that 
r = < (2L — l)*^ (the number of vertices of G cannot be larger than the number of sites in A). 

Equation ()6.25|) cannot be uniformly bounded since, from the recursion relation n(n, r) = 
Il{n — 1, r — 1) + r Il{n — 1, r) (see j^), we have 

n(n, r) > rU{n - 1, r) > r"-''n(r, r) = r""^ , 

which gives a factorial growth 

min(n,|A|) ;r> , , , 

^ ' ' (Cn/2)"/2 if n < |A|/2 



J2 n(n,r)C^> 

r=l 



(C|A|/2)""l^l/2 1^1/2 



after replacing the sum by the term with r = min(n, |A|)/2. 

A sharper upper bound for (|6.23|) may be assumed if one think of G[g] as being given by 



C[G] = sup 



n '^b,b' 



(bb')eE{G') 



(l + dist(6i,M)''^ . (6.26) 



As one varies the partition P of {1, . . . , n}, the graph G, and the path F over it, varies accordingly 
and the decay of ^h^b' in this formula can be useful. We propose that C[g] = Cn,r depends on the 
number of vertices r = |V^| and edges n = \E\ as follows. 

Conjecture 6.6 Let n(?T,,r) = C„^r n(r;,, r). There exist a geometric constant G < 00 such that 

n(n, r) < n(n - 1, r - 1) + C Yi{n - 1, r) (6.27) 
holds for n,r E N, n > r with n(r, r) = G^ . 
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Note that n(n, r) satisfies ()6.27|) with C replaced by r. Assuming ()6.27|) and using that 11(1, 1) 
C and n(fc, I) = if k < I, we have 



^n(n,r) < ^n(n-l,r-l) + C ^n(n-l, 

r=l r=2 r=l 

n~l 

= (l + C)5^n(n-l,r) 

r=l 



(6.28) 



< C(l + C)", 
which leads ()6.25|) to be bounded by 

min(n,|A|) 



Kl<J2{W E n(n,r)<^, (6.29) 



n>l r=l 



where 5' = 4(1 + C)S. 

This concludes the preliminaries and we are now ready to prove Theorem 12. 51 We observe at this 
point that no restrictions about the random variables a^'s has been made beside ()6.16|) and ()6.24|1 
with S small enough. Has Conjecture 16.61 been proved one could work along similar expansions to 
show that (-Ai^uj)"^/^ (-A^^^) {-Al^w)'-^^^ converges to [E (/ - -Dl.u,)""^] with probability 1. 

Proof of the upper bound of \2.1!A Let us recall some facts about the matrix — A^ By equation 
()2.4p it is a positive definite matrix and its square root is well defined. We also have E (— A^ „,) = 
—Al;w = —wAl and, by Lemma i(— A/, „,)~^i^/L^ converges with probability 1 to (— 9^(k))~^ 
exactly as i(— V/, ■ L)~^iy L"^ does. 

In view of this, we can apply Schwarz inequality to the following identity^: 

/ = (E (-A^,^)i/2 (-Ai,^)-i/2) ' < E (- A^,^) E (-Ai,.)-i 

in order to get 

(E(-Ai,^)-i)"' < -WAi, (6.30) 
which implies n <w and concludes our assertion. □ 
Proof of the lower hound of \2.1'J[. From equations (|6.1|) . (|6.3|) and (|6.14p . we have 

(E(-Ai,^)-i)-' = -wAL-w{-ALY''QL{-ALf'^ (6.31) 

= -wV*RlV, 

where Rl '■ M'^^ — > is a matrix whose elements, in view of ()6.2H1 and ()6.29|1 . are bounded by 

with Kl < 8'C I (1 — 5'). Note from equations ()6.3()|1 and ()6.31|) that 9l is a positive matrix. 

^In the following, for any two matrices A and B, A < B means (u, Au) < {u, B u) for all vectors u. 
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Using the isometry operators ()3.14j) and ()3.15p and the fact that iH is the identity matrix 
in Ml^l, we have i (E (-A^,^)-!)"^ it = w {L (iRii^) {L iVi^) with LiV*i^ and LiVi^ 

converging in Lq (V) to the operator d = {d/d^i, . . . , d/d^d)- In addition, we claim that the kernel 
of i (/ — Rl) converges in distribution, as L — > cxd, to the delta function S(r],^) times a. d x d 
matrix g, since the matrix elements of Rl decay faster than l/dist(6, fc')*^ as dist(6, 6') —>■ oo. For 
this, note that i (/ - Rl) = ^ {L"^) if ^ ^ and = C (L"') ii r] = i. 

Whether q is a, diagonal matrix cannot be decided by our estimates. The results from this 
section leads to URli^) . • (?7, — ^ (^j-i + Qij) ^{v^O ^^d this implies 

K = w {1 — g) 

where 1 is the d x d identity and p is a positive matrix satisfying g < 5'C / (1 — 5'). 

□ 
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